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ABSTRACT 

This paper outlines a simple approach to evaluate the atmospheric composition of hot rocky planets 
by assuming different types of planetary composition and using corresponding model calculations. To 
explore hot atmospheres above 1000 K, we model the vaporization of silicate magma and estimate 
the range of atmospheric compositions according to the planet's radius and semi-major axis for the 
Kepler February 2011 data release. Our results show 5 atmospheric types for hot, rocky super-Earth 
atmospheres, strongly dependent on the initial composition and the planet's distance to the star. We 
provide a simple set of parameters that can be used to evaluate atmospheric compositions for current 
and future candidates provided by the Kepler mission and other searches. 

Subject headings: Astrobiology - Atmospheric effects - Planets and satellites: atmospheres - Planets 
and satellites: composition 
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1. INTRODUCTION 

In recent years the search for extrasolar planets has 
resulted in the discovery of the first confirmed rocky 
planets with mi nimum ma sses below 10 (hereafter 
super-Earths) (jUdrv et al. ll2007l : lLeger et al. 112009 ) and 
radi i below 2 Rffi, consi s tent with rocky plan et mod- 
els (jBorucki et al. 1 120111 : iBatalha et 1171120111 ). Many 
of the detected super-Earths are located extremely close 
to their stars, at distances much less than Mercury to the 
Sun. These close-in hot super-Earths can reach temper- 
atures high enough to melt the silicates present on their 
surfaces by irradiation from their host star alone. This 
causes outgassing and the subsequent formation of the 
planet's atmosphere. 

The formation of a rocky planet's atmosphere is a com- 
plex process built up from initial capture and degassing 
of material of me teoritic composition . Our mod el is 
based on work by iSchaefer fc Feglev I (j2004 l2009f ) and 
assumes, as a first order approximation, that at temper- 
atures Tp >1000 K, the material present on the surface of 
the planet vaporizes and forms an outgassed atmosphere. 
In this process, the atmosphere is strongly dependent 
on the planetary com p osition ( Schaefer fc Fegley 20091; 
[Valencia et aTl l2010f ). lSchaeferfc"Feglev I (j2009f ) mod- 
eled the specific case of CoRoT-7b and showed that these 
close-in planets may have a silicate atmosphere as a re- 
sult of outgassing of magma on the surface. 

|miguel@mpia-hd.mpg.de| 



The Kepl er mission recently an nounced 1235 planetary 
candidates (jBorucki et al. I [201 If ), of which 615 are con- 
sistent with models of rocky interior {R < 2zb 0.5 Rq, 
accounting for the 25% error in radii in the Kepler data 
set). 193 of these Kepler objects of interest (KOI) have 
potentially extremely high temperatures on their surfaces 
{Tp >1000 K, when assuming a 0.01 albedo), due to the 
intense stellar irradiation to which they are exposed, be- 
cause of their proximity to the central star. 

In this work, we explore the composition of initial 
planetary atmospheres of hot super-Earths in the Kepler 
planet candidate sample, according to their semi major 
axis. Our approach can be applied to current and future 
candidates provided by the Kepler mission and other ex- 
oplanet searches. 

2. OUTGASSING MODEL 



We use the MAGMA c ode (iFeglev fc Cameron I [19871 : 
ISchaefe r fc Feg levl l2004 l which calculates the equilib- 
rium between the melt and vapor i n a ma^ma exposed 
at tem perature s higher than 1000 K (jHastie et al. IIT9821 
IT981 iHastie Bonneh I IT9851 [19861 ). for Al, Ca, Fe, K, 
Mg, Na, O, Si, Ti and their compounds. 

We assume that the initial composition of a hot super- 
Earth is similar to the composition of komatiites, which 
are ultramafic lavas that were erupted during the Ar- 
chaean on Earth (3.8 to 2.5 billion years ago), when 
Earth had a higher surface tempe rature. This choic e 
is also supported by the work by Ruble et al. I (|2Q04l ). 
who showed, for the case of the Earth and Mars, that a 



Semimajor axis (AU) 

Fig. 1. — Effective stellar temperature vs semi-major axis of all 
Kepler planetary candidates, in different point sizes, according to 
their radii. Planet candidates with Tp >1000 K and Rp < 2.5 Rq, 
are shown in red. The cooler planets (in the same radius range) 
are shown in black. The planet candidates with Rp >2.5 are 
shown in grey. 

more massive planet has a lower composition of FeO in 
the mantle than a less massive one, a consequence of the 
solubility of oxygen in liquid iron-alloy, which increases 
with temperature. Table [1] shows komatiite and bulk sil- 
icate Earth (BSE) compositions (see section I3.1.ip for 
comparison. 

3. ATMOSPHERIC COMPOSITION OF KEPLER 
CANDIDATES 

The planet's equilibrium temperature due to irradia- 
tion by a central star of radius and effective temper- 
ature Teff^^^ is given by equation (TJ 
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where a is the planet's semi-major axis and A 
its bond albedo. We ad opt an albedo of 0.01 
(|Miller-Ricci fc Fortnev~ll2010f ). Different values for the 
albedo can change our results slightly (see discussion). 
Here we set Tp equal to the planet's initial surface tem- 
perature, to explore the resulting atmospheric composi- 
tion. 

Figure [T] shows the stellar effective temperature ver- 
sus planet semi- major axis for all the Kepler candidates, 
where different point sizes indicate different planet ra- 
dius. Using equation [TJ we separate the planets accord- 
ing to their radius and surface temperature. The planet 
candidates with derived Tp > 1000 K and Rp < 2.5 Re, 
are shown in red and are the sample we focus on here. 
The cooler planets (in the same radius range) are shown 
in black. The planet candidates with Rp >2.5 are 
shown in grey. 

3.1. Simulation Results 

We perform outgassing simulations for temperatures 
between 1000 and 3500 K, and apply our results to 
planet candidates in the February 2011 Kepler release 
with Rp < 2.5 Re and Tp > 1000 K. Table El lists their 
semimajor axis, radii, surface temperatures (obtained 
adopting an albedo of 0.01), stellar effective tempera- 
tures and stellar radius. 



Fig. 2. — ^e//,* vs. a for the Kepler candidates analyzed in this 
work (upper panel) . KOIs with radius between 2-2.5 (large 
points), radius between 1.5-2 (medium points) and Rp <1.5 Rq 
(small points). Lines of constant temperatures are shown. Planet 
surface temperatures vs. partial pressures of the gases vaporized 
from a komatiite magma (lower panel). Temperatures at which the 
dominant gases change are indicated. 



Following iSchaefer fc Fegleyl ()2009f ). we explore the 
partial pressures of all the gases vaporized at differ- 
ent temperatures, assuming atmospheres free of volatile 
gases such as H, C, N, S (see discussion). We model the 
initial planetary atmosphere characteristics depending 
on the radius and semi-major axis of the Kepler candi- 



dates. Figure [2( a) | shows the planet candidates that orbit 
their star witn short periods and accordingly high tem- 
peratures. The candidates are shown in different point 
sizes according to their radius: 2-2.5 Re (large points), 
1.5-2 Re (medium points) and candidates with Rp <1.5 
Re (small points). The lines indicate constant planet 
surface temperatures. The b lue, violet, light blue, grey 
and black line in figure |2(a)| indicate Tp of 2096, 2460, 
2735, 2974 and 3170 K, respectively. These temperatures 
indicate significant chan ges in the partial pressures of the 
dominant gases. Figure [2(b)] shows the partial pressures 
of the gases as a function of the planet surface tempera- 
ture. 



Planets located below the blue line in figure 2(a) have 
calculate d sur face temperatures Tp <2096 K. As seen 
in figure 2(b)[ these planets are characterized by par- 
tial pressures dominated by Na, O2, O and Fe (from the 
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TABLE 1 

Compositions adopted for the magma 



Oxide Si02 MgO AI2O3 Ti02 Fe203 FeO CaO Na20 K2O 



Komatiite (%) 47.10 29.60 4.04 0.24 12.80 0.0 5.44 0.46 0.09 

Bulk Silicate Earth (%) 45.97 36.66 4.77 0.18 0.0 8.24 3.78 0.35 0.04 



TABLE 2 

Planetary and stellar parameters for Kepler candidates 

WITH RADII LESS THAN 2.5 R© AND TEMPERATURES HIGHER THAN 

1000 K. 



KOI 


a (AU) 




Tp (K) - 


Teff,. (K) 


(Ro) 


69.01 


0.056 


1.6 


1129 


5480 


1.03 


70.02 


0.046 


1.6 


1001 


5342 


0.7 


72.01 


0.018 


1.3 


1946 


5491 


0.98 


85.02 


0.035 


1.7 


1987 


6006 


1.66 


85.03 


0.084 


2. 


1282 


6006 


1.66 


107.01 


0.075 


2.1 


1025 


5816 


1.01 


112.02 


0.048 


1.7 


UU 


5839 


1.22 


115.02 


0.076 


2.2 


1251 


6202 


1.34 


117.02 


0.058 


1.3 


1142 


5725 


1 


117.03 


0.043 


1.3 


1326 


5725 


1 


123.01 


0.071 


2.3 


1188 


5897 


1.25 


124.01 


0.111 


2.3 


1006 


6076 


1.32 


137.03 


0.046 


2.3 


1335 


5289 


1.27 


139.02 


0.045 


1.2 


1272 


5921 


0.9 



Note. — Table [2] is published in its entirety in the electronic 
edition. A portion is shown here for guidance regarding its form 
and content. 

^ Calculated adopting an albedo of 0.01 

highest to the lowest partial pressures). In the same way, 
planets located between the blue and the violet line have 
temperatures 2096 < Tp < 2460 K. In this case, the 
partial pressure of SiO becomes higher than Fe. With 
surface temperatures between 2460 and 2735, the planets 
are located between the violet and the light-blue line. For 
these surface temperatures, the partial pressure of SiO is 
higher than O. When the surface temperature is between 
2735 < Tp < 2974 K (between the light-blue and the 
grey line), the partial pressure of SiO is higher than O2 
and for 2974 < Tp < 3170 K (between the grey and the 
black line) SiO becomes more abundant than monatomic 
Na. Finally a planet with > 3170 K is located above 
the black line and is characterized by the highest partial 
pressure of SiO, followed by O2, monatomic O and Na. 

3.1.1. Different types of Atmospheres 

In order to link the partial pressures to the result- 
ing atmospheric composition, we calculate the column 
density (cr^) of each g as as a function of the gas' 
partial pressure, P^, for planetary masses between 1 
and 10 Mffi, corresponding to radii between ^ 1-2.5 
Rffi ("Valencia. Sasselo v fc O'Connell II2007I : iSeager et al I 
[2007 : Sotin et al. 20071 ) using equation El 



(2) 



where Na is Avogadro's number, /i^ is the molecular 
weight of the species and g is the gravitational surface 
acceleration of each planet. We find that the composi- 
tions of the atmospheres do not show significant changes 
when considering different planetary masses. This means 



that the surface temperature and therefore the distance 
to the host star is the dominating factor in determining 
the composition of the atmosphere, assuming a certain 
planetary composition. 

We adopt a komatiite composition for hot super-Earths 
and explore the effects of initial composition on the re- 
sults by comparing with a Bulk Silicate Earth (B.S.E) 
initial c ompo sition in Figure [3l Figure |3(a) (komati- 
ite) and 3(b)| (BSE composition) show the stellar effec- 
tive temperature vs. the planet's distance to the star. 
Lines show the temperatures that separate the different 
classes of planetary atmospheres. Areas highlighted in 
gray scale indicate the regions of different types of atmo- 
spheres. The planets in each one of these areas are char- 
acterized by di fferen t t ypes o f atmospheres as indicated 
below. Figures 3 (c) | and |3 (d) | show the column density of 
all the gases vaporized vs. the planet surface tempera- 
ture for a 10 M0 planet with komatiite and BSE as the 
initial composition, respectively. Since we find similar 
results for planets with different masses, the column den- 
sities shown in Figures [3(^ and [3(d)] are representative 
of planets with masses between 1-10 Mq. The dominant 
gases are Fe, Mg, Na, O, O2 and SiO. 

Our results divide hot rocky planet atmospheres into 5 
classes. Temperatures that limit each class change when 
considering different magma compositions. For komatiite 
composition: 

Atmospheres Type I {Tp < 2033 K), characterized 
by the presence of monatomic Na, O2, monatomic O and 
monatomic Fe in order of abundance. 

Atmospheres Type II (2033 < Tp < 2588 K), 
SiO becomes more abundant than monatomic Fe. The 
atmosphere is characterized by monatomic Na, O2, 
monatomic O, SiO, monatomic Fe and monatomic Mg 
from the most to the least abundant. 

Atmospheres Type III (2588 < Tp < 2890 K), 
monatomic Mg becomes more abundant than monatomic 
Fe. The gases with the highest column densities are 
monatomic Na, O2, monatomic O, SiO, monatomic Mg 
and monatomic Fe. Note that the column densities of 
O2, monatomic O and SiO are almost equal in this tem- 
perature range. 

Atmospheres Type VI (2890 < Tp < 3168 K), SiO 
becomes more abundant than O2 and the atmosphere 
becomes silicate and monatomic Na dominated. 

Atmospheres Type V {Tp > 3168 K), dominated by 
silicate oxide, followed closely by monatomic O, O2 and 
monatomic Na. 

Of the 193 Kep ler pl anetary candidates studied in this 
work (see Figure 3(a)), 189 are characterized by Type I 
atmospheres when adopting komatiite as the initial com- 
position. Three planets are characterized by Ty pe II (be- 
tween the black and the dark-red line in figure p^aj), no 



planets by Type III (between the dark-red and red lines), 
only one by Type IV (located between the red and yel- 
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Fig. 3. — Stellar effective temperature vs. semi-major axis of Kepler planetary candidates (upper panels, see figure [2(a)| for detailed 
explanation) and column densities of the atmospheric gases vs. planet surface temperature (lower panels ). Th e colunin densities were 
calculated for a planet of 10 Mq, but are r epresentativ e of planets with masses between 1-10 Mm (see sec tion ts.l.ll) . We adopt two different 
planetary compositions: komatiite (figures |3(a)||3(c)| ) and bulk silicate of the Earth (figures [3(b)l |3(d)| ). 



low lines) and none by Type V atmospheres (above the 
yellow line) , if we only consider the irradiation from the 
star for heating (see discussion). 

When adopting the BSE composition, the temperature 
limits between the 5 classes change slightly to 2070, 2120, 
2996 and 3278 K, respectively. For a BSE composition, 
190 planetary candidates are characterized by Typ e I at - 
mospheres (located below the black line in figure 3(b)), 
only 1 planet by Type II atmosphere (between the black 
and the dark-red lines), 1 by Type III (between the dark- 
red and red lines), 1 by Type IV (located between the red 
and yellow lines) and none by Type V (above the yellow 
line). As a result, note that, for this planetary sam- 
ple, the initial composition is not important when the 
planet's temperature is high for atmospheres of Type IV 
and V. Nevertheless, for lower temperatures (Types I, II 
and III) it becomes an important factor, since a small 
difference in temperature can change the dominant gases 
in the planetary atmosphere. 

4. DISCUSSION: THE SCOPE OF OUR WORK 

The composition of the atmosphere depends strongly 
on the composition of the planetary interior. We show 
the changes in the characterization when adopting two 



known compositions taken from our Solar system, but 
substantially different compositions are possible. 

Following Schaefer & Fegley (2009), we consider at- 
mospheres produced by the outgassing of a volatile- 
free magma. At these extreme temperatures, elements 
such as H, C, N and S should escape from the at- 
mosphere in short timescales and consequently are not 
present in the atmosphere. This assumption is sup- 
ported by observations of volatile loss in rocky bodies 
(■Spencer & Schneider 1996; Fegley 2004) and studies in 
h ot-Jupiters (Bara ffe et al. 2005) and hot super-Earths 
(|Leitzingers et al. 11201 If ). If the volatiles were not lost, 
the atmosphere would be dominated by volatiles. In 
this paper we only concentrate on close-in, hot plan- 
etary candidates, consistent with volatile free atmo- 
spheres. Cooler planets that could potentially provide 
habita ble conditions ar e not discussed here (see e.g. 
Kalte negger fc Sasselovl (j2011[ )). 

In this work we have adopted a bond albedo for the 
planets of A =0.01. This a ssumpti on is suppor t ed by 
obser vations ([Rogers et al. 2009; iRowe et al. I l2006l 
2008 ) and modehng (Sudarskv. B urrows fc Pintol 
20001: [Burrows Ibgui & HubenyT 120081 : 



Miller- Ricci fc Fortnev 1 1201 0[ ) of hot exoplanets, which 
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find very low albedo. Although there are currently no 
indications that suggest higher albedos for hot rocky 
planets, we studied the effect of a higher albedo on our 
results. When changing the albedo from 0.01 to e.g. 
0.2, the size of the sample is affected. For an albedo of 
0.2, only 146 planetary candidates with radius less than 
2.5 R® have Tp > 1000 K. Of these 146 planets, 144 are 
characterized by atmospheres Type I, 1 by Type II, 1 
by Type III, by Type IV and by Type V (adopting 
komatiite composition) . 

We assume here that the planets are only heated by the 
irradiation of the star. Other heat sources like short-lived 
radioisotopes, impacts, core formation and tidal effects 
can in crease the surface temperature. iHenning et al. I 
([2009) showed that the insolation is the most important 
heating source for planets with periods larger than ~ 2 
days and for planets with shorter periods and non ec- 
centric orbits. For very eccentric planets (e >0.1) which 
orbit their star with periods less than ~ 2 days, tidal 
heating dominates over insolation. We explore the ef- 
fects of tidal heating in the Keple r planetary candidates 
with P < 2 days, using eq. [3] (jHenning et al. 1 120091 : 
iMurrav fc Dermott 1120051 



' p,tidal-\-ins 



nRlil - A) 



47ra^ 



21 fe: 



GRlne' 



Mt 



(3) 

where, k2 is the second-order Love number of the planet. 
For a super-Earth we adopt 0.3 (Henning personal com- 
munication), Q is the quality factor of the planet, n its 
mean motion and e the eccentricity, which is unknown 
for the Kepler candidates. We adopt different values for 
e and Q. Most of the candidates that might be affected 
are located very close to the star (a less than ~ 0.03 
AU), where we expect almost circular orbits. We adopt 
e =0.001 and e =0.01 and two different values for Q: 200 
and 20000, which represent planets with limited partial 
melting or significant parti al melting in the mantle, re- 
spectively (IHenning et al. '^09), to explore this effect. 
In the four cases analyzed, we find that there is no signif- 
icant change in the surface temperatures of these planets 
when considering heating by insolation plus tidal effects 
compared to heating by insolation only. The role of tides 



in these planets, if occurring, will be to increase the con- 
vective vigor of the mantle, causing a larger fraction of 
total mantle volume to be accessible to surface degassing, 
but will not change the surface temperature significantly. 

5. CONCLUSIONS 

In this work we explore the different atmospheric com- 
position of hot super-Earth candidates from the February 
iTep/er database according to their radii, semi-major axis 
and stellar effective temperature. 615 Kepler planet can- 
didates have radii less than 2.5 Re, of which 193 have 
calculated surface temperatures larger than 1000 K (as- 
suming an albedo of 0.01), and are the focus of this work. 
We model the initial planetary atmospheric compositions 
for komatiite and bulk silicate Earth magma composi- 
tions. The planet surface temperatures were calculated 
assuming insolation only. We discuss the consequences of 
adopting different sources of heating and find that tidal 
heating is not important when adopting low eccentrici- 
ties, consistent with the Kepler sample. 

The results indicate 5 types of atmospheres, based on 
the abundance of the dominant gases (Fe, Mg, Na, O, 
O2 and SiO). When adopting komatiite composition, 189 
of 193 Kepler candidates present a Type I model atmo- 
sphere dominated by monatomic Na, O2, monatomic O 
and monatomic Fe. Three planetary candidates are char- 
acterized by Type II atmospheres, where SiO becomes 
more abundant than monatomic Fe, no candidates with 
Type III atmospheres, where monatomic Mg is more 
abundant than monatomic Fe, and only one planet is 
characterized by Type IV, with a large abundance of SiO. 
Finally, there is no Kepler planetary candidate with an 
atmosphere of Type V, which is completely dominated 
by SiO. 

This atmospheric characterization is almost indepen- 
dent of the planetary mass between 1-10 Mq and low 
eccentricities. The initial magma composition is impor- 
tant for the resulting atmospheric gases, especially for 
lower surface temperatures. We find that the distance 
from the host star is the dominant factor to characterize 
a rocky planet's atmosphere. 

We are grateful to Dimitar Sasselov and Wade Henning 
for stimulating discussions. 
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